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Abstract 

The notion of fractional dynamics is related to equations of motion with one or 
a few terms with derivatives of a fractional order. This type of equation appears 
in the description of chaotic dynamics, wave propagation in fractal media, and field 
theory. For the fractional linear oscillator the physical meaning of the derivative 
of order a < 2 is dissipation. In systems with many spacially coupled elements 
(oscillators) the fractional derivative, along the space coordinate, corresponds to a 
long range interaction. We discuss a method of constructing a solution using an 
expansion in e = n — a with small e and positive integer n. The method is applied 
to the fractional linear and nonlinear oscillators and to fractional Ginzburg-Landau 
or parabolic equations. 

PACS: 45.10.Hj; 45.05.+x; 45.50.-j 
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1 Introduction 

It became clear in the last decade that many physical processes can be adequately described 
by equations that consist of derivatives of fractional order. In a fairly short period of time, 
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the list of such apphcations is long and the areas of applications are broad. Even in a concise 
form, the applications include material sciences [1, 2], chaotic dynamics [3], quantum theory 
[4, 5, 6, 7], physical kinetics [3, 8, 9, 10], fluids and plasma physics [11, 12, 13, 14], and many 
others physical topics related to anomalous diffusion, wave propagation [15], and transport 
theory (see more in reviews [3, 16]). Since the fractional calculus has a fairly long history, the 
approaches for solutions of corresponding equations are also enormously rich. Let us mention 
some of these approaches that are related more specifically to this paper: probabilistic basis and 
interpretation of the fractional kinetics [17, 18, 9, 19], dissipative interpretation of the fractional 
derivative [20, 21, 22, 23], Green's function method [24, 25], etc. 

This paper is motivated by a lack of methods that permit an explicit analysis of the equations 
with fractional derivatives that describe some important physical processes. Particularly, the 
difficulties are due to the absence of characteristic scales (sometimes such processes are called 
multiscaling) and the difficulty of direct estimates. In typical physical situations the "level of 
fractality" is low, i.e. the order of fractional derivatives a of the corresponding terms of the 
dynamical equations slightly deviates from an integer value n (in the considered cases n = 1 or 
2) . This recalls a possibility to use an expansion over the small parameter e = n — a that we 
call an e-expansion. We develop a construction that is applied to linear and nonlinear fractional 
oscillators. Particularly, for the linear fractional oscillator (LFO) the obtained expansion can 
be compared to the exact one. As a more complicated example, we consider a solution of the 
fractional Ginzburg-Landau equation (FGL) introduced in [26, 27] (see also [28]). 

The basic description of the £-expansion is given in Sec. 2 for different fractional derivatives. 
In Sec. 3, we consider examples of the £-expansion application to linear and nonlinear oscillators. 
In Sec. 4, we consider large t asymptotics {t oo). Application of the e-expansion to the 
Landau-Ginzburg or to the parabolic nonlinear equation is given in Sec. 5. Some technical 
details are the Appendices 1, 2. 
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2 Description of the ^-expansion 



2.1 Caputo fractional derivative of order ol = 2 — e 

The fractional derivative has different definitions [29, 30], and exploiting any of them depends 
on the kind of the problems, initial (boundary) conditions, and the specifics of the considered 
physical processes. The classical definition is the so-called Riemann-Liouville derivative [29, 30, 
31] 



r(n - a) ax" J a {x - ' 

* ''•^^ V{n-a)dx^]^ (^_^)a-n+i' 
where n — 1 < a < n. Due to reasons, concerning the initial conditions, it is more convenient to 
use the Caputo fractional derivatives [2, 32, 33]. Its main advantage is that the initial conditions 
take the same form as for integer-order differential equations. 

The left Caputo fractional derivative [2, 32, 33] is defined by the equation 

I ft fW(r) 
^ ^ ° * ^ ^ T{n-a) Jo {t- r)°+i-" ' ^ ' 

where n — 1 < a < n, and f^^'^r) = (f^ f (r) / dT'\ For n = 2, 
This presentation is not convenient in the limit e — 0, since 

^ =8 + 0(8% (4) 



r(£) 

and we present (2) and (3) in the form [25]: 
/(")(0)f 



DVit) = ^ + TT. -TT / /^"+'Hr)(i - r)--"dr, n-l<a<n, (5) 

L in — a + 1) i (n — a + Ij Jo 



r{n -a + 

or for n = 2, 



•'^ ^ r(i + £) r(i + e) io 

Let us consider first the case 

£ i < 1. (7) 
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We can use the expansion 



i (t - Tf = — e^Ht-r) = 



r(l + £)' ' r(l + e) 

= 1 + £ (7 + ln(i - r)) + £2 (^^ \n\t - r) + 7ln(i " ^) + y " ^) + -> (§) 
where t < t, and 7 = 0.5772156649... is a constant. As the result, we obtain 

D'-'f{t) = f^'Ht) + e (^/(2)(0) ln(t) + 7/(2) (t) + /(3)(^) in(t - r)dr) + ... (9) 

that expresses the fractional derivative in a form of the perturbation to the second derivative, 
when £ ^ 1: 

D^-^f{t) = f^'\t)+EDlfit) + .... (10) 

where Df is defined as: 

Dffit) = /(2)(0) ln{t) + jf(''\t) + f f(^\T) \n{t - T)dT. (11) 

■/o 

Note that, the limit £ — >^ in equation (9) gives the correct expansion 

Another useful comments is that the 3d derivative J^^\t) should exist in order to use the 
correction of order e in (9). 

As a simple example of application of formula (9), consider the a-derivative of and t^: 

D'^-H'^ = 2 + 2£7 + ... 

D'^-H^ = Qt + £[6(7 - l)t + 6t ln(t)] + ... (12) 
Prom the well-known relation 



we get the exact results 



na^H _ r(/3 + 1) 



r(i + £)'' 

r(2 + £) ^ ^ 



Expansions of (14) for £ t ^ 1 coincide with (12). 



2.2 Caputo fractional derivative of order a = 1 — e. 



To consider the limit a —> 1, let us put n = 1 in (5): 

Similarly to n = 2, consider the case e t ^ 1, and use the expansion (8). It gives 

D'-'f{t) = f^'\t)+eDlfit) + ..., 

where Dlf{t) is 

Dlfit) = /(^)(0) ln(t) + 7/(1) (0 + /* /(2)(r) ln(t - r)dr. 

Jo 

Note that, if we consider the limit s ^ 0, we get 

\im D'-^ f it) = f^'\t). 

2.3 Riesz fractional derivatives of order a = 2 — e 

The Riesz fractional derivative of order a is defined by 



d\x\°' ^ ' 2cos(7ra/2) 
where a 7^ 1,3,5..., and are Riemann-Liouville fractional derivatives 



r(n — a) ax'^ A 
Substitution of Eqs. (20) into Eq. (19) gives 



r(ra - a) dx'' J-00 [x - z)"-"+i ' 

l)n Qn .00 /(2)d^ 



_ ^\a— n+1 ■ 



(Z-X) 



d\x\<^ ^ ' 2cos{Tra/2)T{n - a) dx"- \J-oo {x - z)'^-"'+^ Jx (z - a;)«-"+i 
For the order a = 2 — e, 

2 cos(7r£/2)r(£) Sx^ Vi-oo (x - z)^'' ^ (^ - x)i-^ 

Using the condition /(itoo) = 0, we get 

r = r mix-zr-'dz = l r f('\z)ix-zrdz, 

J —00 yX J —00 £ J —00 



f{z)dz _ r j^^^^^_^y-i^^^_l r f{^)^z){z-xfdz. (23) 



{z - xy-^ Jx e 
These relations and equation £r(£) = r(l + e) allows us to present (21) in the form 

fja 1 fl2 / i-x roc 

fi^) = / f^'K^K^ - ^fdz - / f^'\z){z - xfdz ) . (24) 



d\x\'^'^ ' 2cos(7r£/2)r(l + £)ax2 
For the case e a; ^ 1, we use 

{x — z) 



oo 



cos(7re/2)r(l + e] 

to obtain 



= l + £[7 + ln(a;-z)] + ... (25) 



"^'Z^fi^) = f'-'H^) + eDl^fix) + (26) 



d\x\'^-^' 
where 

Dmfix) = lf^^\x) + [/(3)(x - y) - /(3)(x + y)\ \n{y)dy. (27) 

A general comments to this section is that all results (11), (17) and (27) are valid under 
restriction £t ^ 1. The expansion of Riemann-Liouville derivative is considered in Appendix 1. 
For the asymptotic t — ^ oo, and 

£t > 1 (28) 
another expansion is necessary. The corresponding results will be considered in section 4. 



3 Linear and nonlinear fractional oscillator 
3.1 £-expansion for linear oscillator 

A linear fractional oscillator (LFO) is defined by the equation 

r>°a;(i)+a;^a;(i) = 0, (29) 

where uj is dimensionless "frequency", a = 2 — £, 0<£<Sl, and D"^'^ is Caputo fractional 
derivative that allows us to use the usual initial conditions [25]. LFO is an object of numerous 
investigations [20, 21, 22, 23, 34, 35, 36, 37, 38] because of different applications. 
The exact solution of Eq. (29) for 1 < a < 2 is [20, 21]: 

x{t) = x{0)Ea,i{-ioH'') + tx'{0)Ea,2{-ioH''), (30) 
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where 

is the generahzed two-parameter Mittag-LefHer function, and 

oo 

EJz) = Eaiiz) = y — — — 

i^r(aA; + l) 

is one-parameter, or simply, Mittag-Leffler function [39, 40]. 

Existence of an exact solution of (29) in the form (30) permits a comparison with 
expansion 

x{t) = xo{t) + exi{t) + ... 
constructed below. The equation for xo{t) and the corresponding initial conditions are 

x'o'(i) +u'^xo{t) = 0, xo(0) = x(0), 4(0) = x'{0), 

with the solution 

Xo{t) = x{0) cos{u!t) + {x'{0)/u!) sin(a;t). 
The equation for xi (t) and the initial conditions are 

x'l{t) + uj'^xi{t) + Djxoit) = 0, xi(0) = 0, x'i(O) = 0, 

where Z)fxo(t) is defined by (11): 

Dfxoit) = x'^{0) ln{t) + jx'^it) + / x^o\t) ln{t - T)dT. 

Jo 

Substitution of Eq. (35) into Eq. (37) gives 

DiXo{t) = -u'^x{0) \n{t) - 7a;2[x(0) cos(a;t) + {x'{0)/u) sin(a;i)] + 

+u^x{0) f sm{ujT)\n{t - T)dT - uj'^x'{0) f cos(u;t) ln(i - t)o?t. 
Jo Jo 

The first and second integrals can be calculated, and we have the relations 

rt 



/ sin(a;r) ln{t — T)dT = 
Jo 

= uj~^ ^ln(t) + cos(a;i) [ln(a;) + 7 - Ci(a;t)] + sin(a;t) 



7 + 1 ■ 

Si(wt) 

2 4 ^ ' 



and 



/ cos(wr) ln(i — r)dT 
Jo 



= a; ^ ( cos{ujt) 



7+1 TT 



- SiM) 



+ sin(a;i) [- ln(a;) - 7 + Ci(a;t)]^ . 
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Here Ci(t) and Si(t) are sine and cosine integral functions [42] respectively: 

, /■* sin(x) , ... ,„ X f°° cos(x) , ... 

Si(t) = / — ^dx = si(t) + 7r/2, Ci(t) = - / ^dx = -ci(t). 

Jo X. Jt X. 
As the result, Eq. (36) transforms into 

Xi{t) + u'^xi{t) + u}^x{0) cos{ut) + uvo sin(a;i) (ln(u;) — Ci(u;t)) + 



(40) 



+ 



a;^x(0) sin(u;t) — lovq cos^cvt) 



^-^-Si(a;t))=0. 



(41) 



To simplify the calculations consider x(0) = 1, x'{0) = 0, and u = 1. Then xo(t) = cos(t), and 

'7 + 1 TT 



x'({t) + xi{t) - cos(t)Ci(t) - sin(t)Si(t) + ( - sin(t) = 0. 



(42) 



The solution of this equation is 



xi{t) = cos(t) / 

Jo . 



7+1 TT 

^ 4 



[1 — cos^('u)] — sin('u) cos('u)Ci(u) — Si(u) + Si(u) cos^(u) 



+ sin(t) 



— ) sin(ti) cos(w) + cos^(ti)Ci(u) + sin(u) cos(u)Si(u) 



which after integration takes the form 

xi{t) = [sin(i)Ci(t) + cos(t)Si(t)] - ^ + ^ cos(t). 
FinaUy, Eq. (33) provides the solution 

x{t) = cos(t) - |(l - cos(t) + 1 [sin(i)Ci(t) + cos(t)Si(t)]) + ... 



du+ 

(43) 

(44) 



(45) 



The term, that is proportional to e, can be considered as a correction to the solution x{t) = cos(t) 
due to the fractional derivative of order a = 2 — s. 

Let us compare this result with the exact solution for FLO. The decomposition of (30) is 
[20]: 

X{t) = X{0) [fa,o{t) + 9aAt)] + IfaAt) + 9a,l{t)] , (46) 
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where 



TT ^0 



roo r"-!-'' sin(7ra) 



Q-^^ ■ — ^ ^ dr 

j.2a + 2r" cos(7ra) + 1 ' 



5a,feW = ^e*'°'^''/"^ cos[tsin(7r/a)-7rA;/«], (fc = 0,l). (47) 
For the initial conditions a;(0) = 1, and x'(0) = 0: 

x{t) = E^{-r) = fa,o{t) + 9a,o{t)- (48) 

Using a = 2 — £ and the condition e f ^ 1, we get 

fa,o{t) = ~| + 1 sin(i)Ci(t) + t cos(t)Si(t) - t^ cos(t)^ + 

ga,oit) = cos(i) + £ - cos(i) + .... (49) 
Substitution of (49) into (48) gives Eq. (45). 

3.2 Fractional oscillator in phase space 

The variety of definitions of fractional derivatives causes some inconvenience in their applications 
to different physical problems. Particularly, this comment concerns the phase space definition 
for the fractional dynamics, which is not uniquely defined in the general case. As an example, 
consider the fractional Caputo derivative of order a = 2 — e that can be defined as 

where is a fractional integration of order e. Therefore, we can use the relation 

D^-' = J'^D'^ = J'^D^D^ = {J'D^)D^ = D^~'D^ (51) 
to represent the fractional equation 

D'^-'x{t) = F{x{t)) (52) 

as a system of two equations 

D'x{t)=p{t), D'-'p{t) = F{x{t)), (53) 
where first equation has derivative of integer order. 
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In general, we have an inequality 

^l-£l^l-£2 _^ ^2-(£l+£2)^ (54) 

where D^~^ = J^D^, i.e., the composition doesn't exist for the Caputo derivatives 

D^-''D^-'^x{t) = L»2-(^i+^2)^(i) + 1+^2-1^/(0), (55) 

r(ei + 82) 

and the operations and do not commute: 

+£-1 

D^J^x{t) = J^D^x{t) + :p-r^x(0). (56) 
The symmetrized fractional equation 

D^-'/^[D^-'/^x{t)] = F{x{t)), (57) 

is not the same as (52), since 

^l-e/2pl-£/2^/^^^ ^ D^-e^U^ + ^^^'(0). (58) 

r(£) 

Eq. (57) is equivalent to (52) for x'{0) = 0. Consider Eq. (57), and represent it as a system of 
fractional equations 

D^-'/^x{t) = p{t), D^''/'^p{t) = F{x{t)), (59) 
where both equations are fractional equations. Using Eq. (16), we can rewrite (59) as 

x'{t) + {e/2)D\x{t) + ... = p{t), p'{t) + {e/2)D\p{t) + ... = F{x{t)), (60) 

and consider their solution in the form 

x{t) = XQ{t) + exi{t) + p{t) = po{t) + epi{t) + ... (61) 

As the result, we get a system of equations: 

4it)=Po{t), p'oit) = F{xoit)), 

x[{t) + Dlxoit) = pi{t), p[{t) + Dlpoit) = f^^\ a;i(0. (62) 

V / X=XO 

For the initial conditions xo{0) = 1, po(0) = 0, and F{x) = —x, the perturbation terms DIxq 
and -Djpo in (62) are defined as 

rt 



DlxQ{t) = — 7sin(t) — / cos(a;r) ln{t — T)dT, 
Jo 
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Dlpo{t) = - ln(t) - 7cos(t) + / sin(wT) ln{t - T)dT 

Jo 

and leads Eq. (62) to the form 

"7+1 TT 



x[{t) + cos(i)Si(t) - sin(t)Ci(t) - 
p[{t) - cos(t)Ci(t) - sm(t)Si(t) + 



2 4 

7 + 1 tt' 



cos(t) = pi(t), 
sin(t) = — xi(t) 



with solution 



xi{t) = -fsin(t)Ci(t) -tcos(t)Si(t) +cos(t) - 1, 
p^(t) = -sin(t)Ci(t) - cos(t)Si(t) - t cos(t)Ci(t) + t sin(t)Si(t) - sin(2t) - sin(t). 
This solution describes the phase space evolution of FLO for et <^1. 



(63) 



(64) 



(65) 



3.3 Nonlinear fractional oscillator 



Let us consider a fractional nonlinear oscillator that is defined by the equation 



D^-^xit) - x{t) + 2x%t) = 0, 



(66) 



where < e <C 1, and D is a Caputo fractional derivative. For the soliton-type solution 



xo{t) 



1 



sech(t), 



cosh(i) 

the correction of the order e can be found similar to (45): 

x'l{t) + [Qsech^{t) - l]xi{t) + Dlsech{t) = 0, 



(67) 



(68) 



where 



Dlsech{t) = - ln{t) + 7- 



cosh^ {t)-2 /■* sinh^ (r ) [cosh^ (r) - 6] 



cosh-^(t) 



+ 



cosh^ (r) 



ln{t - T)dT. (69) 



The solution x\{t) of equation (69) can be derived numerically. As the result we can get the 
corrections to the solution xo(t). 
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4 Asymptotic behavior for large t 



4.1 Asymptotic representation of fractional derivative 

In the previous sections, we consider the intermediate time asymptotic Umited by the condition 
et <^1. In this section, we consider the opposite condition t ^ oo. 

It is known [32, 25], that the Laplace transform of the Caputo fractional derivative is 

I n— 1 

e"^* [D'^xit)] dt = - J2 s""''"^a;('=)(0), (70) 

k=0 



f 

Jo 



where n — 1 < a < n, and X{s) is a Laplace transform of x{t): 

roo 

X{s) = / e-'^x{t)dt. (71) 
Jo 

Note that formula (70) involves the initial conditions x^''\0) with integer derivatives x^''\ There- 
fore we can put the initial conditions in the usual way. The functions x{t) satisfy the condition 

roo 

/ e-'^\x{t)\ < oo. (72) 
Jo 

For 1 < a < 2, Eq. (70) has the form 

/•oo 

/ e-^* [D'^xit)] dt = s"X(s) - s^-^xiO) - s'^-^x'{0). (73) 
Jo 

Inversion of (73) gives 



1 



D'^x{t) = ^ I e'' s"X(s) - s"-ix(0) - s'^-'^x'iO) ds. (74) 



27ri J Br 

where Br denotes the Bromwich contour, i.e., a line from c — ioo to c + ioo, where Re{s) = c, 
and the contour is taken to the right of all singularities in order to insure condition (72). Closing 
the contour to the right will yield x{t) = for t < 0. For non-integer a the power function s°' is 
uniquely defined as s" = |s|"cxp[i arg(s)], with — tt < arg(s) < vr, that is in the complex s-plane 
cut along the negative real semi-axis. 

The asymptotic expansion {st S> 1) can be formally obtained by the expanding X{s) in 
powers of s, and then inverting term-by-term. The asymptotic for large t corresponds to small 
s. Consider the Taylor series 

oo 

Xis) = Y.-^x(''\0)s\ \s\<l. (75) 
k=o'^- 
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Substitution of (75) into (74) yields 

27rZ J Br 27rZ J Br 

°° 1 ^ r 



k=0 

Using the analytical continuation of 

to the half-plane Re{f3) > 0, we arrive to the equation 

with the leading asymptotic term 

/-2+£ f-l+e 

D-'-^xit) « -^(0) r(_i^^) - ^'(0) ^ £x(0)t-2+- - £x'(0)^-l+^ {et » 1). (79) 

It is seen from (79) that for arbitrary small a;'(0) the first term can be neglected for sufficiently 
large t. 



4.2 Asymptotics for linear fractional oscillator 

Let us consider some applications of the results of the previous section. For the linear fractional 
oscillator 

D'^xit) +u^x{t) = 0, (80) 

with a = 2 — e, we obtain from (79), 

a;(0) x'(0)t-i+= x(0) n,, x'(0) i,^ , 

^ ^ ~ 2 w 1^ 2 rf ^ - -e^t-^^' + e^r^+' , » 1 . 81 

tu^ r(— 1 + e) uj^ 1(e) tu^ o;^ 

This result can be also derived from the asymptotic of the Mittag-Leffler function using the 
exact solution of (80), 

x{t) = x{0)Ea{-coH'') + x'(0)t^a,2(-w^t"). (82) 
Let us use the following integral representation 
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where Ha denotes the Hankel path, a loop which starts from — oo along the lower side of the 
negative real axis, encircles the circular disc |.^| < 1^1^/" in the positive direction, and ends 
at — cx) along the upper side of the negative real axis. By the replacement ^" ^ ^ Eq. (83) 
transforms into [25, 44]: 

Ea{z) = / I (1< a < 2), (84) 

27na J-y(a,S) t.- z 

where 7ra/2 < 5 < min{Tr , wa} . The contour j{a,6) consists of two rays = {arg(^) = 
—S, 1^1 > a} and S+s = {arg(^) = +S, \^\ > a}, and a circular arc Cs = {|C| = 1, — <5 < CLrc{$,) < 
S}. Let us denote the region on the left from 7(a, (5) as G~{a,S). Then [44]: 

oo 

n=l ^ ' 

and S < I arg(z)| < tt. In our case, z = — a;^t", arg(2;) = vr, and 



In a similar way, we get 



We arrive at the asymptotic result (81) that exhibits an algebraic decay for t oo. This 
algebraic decay is the most important effect of the non-integer derivative in the considered 
fractional equations, contrary to the exponential decay of the usual damped-oscillation and 
linear relaxation phenomena. 

4.3 Nonlinear fractional oscillator 

The results of the previous section can be extended to the nonlinear fractional oscillator. Con- 
sider the equation 

D'^-^x{t) -x{t) + {t) = 0, (88) 
and try to find its solution in the form 

x{t) = xo{t) + exi{t) + ... (89) 

For £ = 0, xo{t) satisfies the equation 

x'^{t) - xo{t) + xl{t) = 0. (90) 
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Consider a particular solution 

xo{t) = \/2sech(t), 

and its Laplace transform [41]: 

X{s) = / xo{t)e~'^dt = 

Jo v2 

where "^{z) is the digamma function 

T'{z) d 



1 ^ fs + 3\ 1 /s + 1 



V2 V 4 



For z < 1 



and therefore 



*(2) = -^-7+y2 + 0(22), 



s + 0(s2), (s^O), 



where 



*(n,2;) 



is the n-th polygamma function, and 



-1.831931188. 



Substitution of (79) into Eq. (88) gives 



x^{t) - x{t) - 



V2 



r(i-a) 



^ 0. 



The leading asymptotic term is 



^(t) -4(t) 2r(l - g) 
^^'^ ^ 6r(l - a) + A(i) ' 



where 



A{t) = (^108^2 t"" + 



-12r2(l - a) + 162t 



-2a 



l/2\ 1/3 



Ci + est"" (t 



and 



ci = 2i/3(_3r(l - a))V6^ C2 = 25/6i8(_3r(i _ a))-^/3_ 
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5 Fractional Ginzburg-Landau (FGL) equation 



5.1 Appearance of fractional derivatives in Ginzburg-Landau 
equation 

Let us recall the appearance of the nonhnear parabolic equation [45, 46, 47, 48], and the FGL 
suggested in Ref. [26] (see also [27, 28]). Consider wave propagation in some media and present 
the wave vector k in the form 

k = ko + K = ko + K|| + k;_l, (102) 

where ko is the unperturbed wave vector and subscripts (||,-L) are taken respectively to the 
direction of ko- A symmetric dispersion law u = Lo{k) for k <C fco can be written as 

iu{k) = io{\k\) « u;{ko) + vg (|k| - A;o) + ^< (|k| - kof, (103) 



where 



and 



k=ko 



1 



|k| = |ko + k| = (ko + K||)2 + KlKiko + K|| + -^n±- (105) 
Substitution of (105) into (103) gives 

uj{k) a! wo + VgK\i + + (106) 

where = oj{ko). Expression (106) in the dual space (" momentum representation" ) corresponds 
to the following equation in the coordinate space 

i^_„Z_,„,g_iA,z4A.Z (107) 

with respect to the field Z = Z{t,x,y,z), where x is along ko, and we use the operator corre- 
spondence between the dual space and usual space-time: 

Uik) ^ K|| ^ 

— = - (-1,)^ — = (108) 
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A generalization to the nonlinear case can be carried out similarly to (106) through a nonlinear 
dispersion law dependence on the wave amplitude: 



u = uj{k, \Zf) Lo{k, 0) + b\Zf = Lo{\k\) + b\Zf (109) 

with some constant b = dLo{k, \Z\'^) / d\Z\'^ at \Z\ = 0. In analogy to (107), we obtain from (106), 
and (108): 

i— = u{k,)Z - iv,— - ^A^Z - ^A||Z + b\Z\^Z. (110) 

This equation is known as the nonlinear parabolic equation [45, 46, 47, 48]. The change of 

variables from (i, re, y, z) to (t, x — Vgt, y, z) gives 

f)7 11 v' 

= ^^^^ + f ^11^ - '^(M^ - b\Z?Z (111) 

that is also known as the nonlinear Schrodinger (NLS) equation. 

Wave propagation in a media with fractal properties can be easily generalized by rewriting 
the dispersion law (106), (109) in the following way [26]: 



(fc, |Z|2) = a;(A:o, 0) + v^k. + 51 («!)"/' + 92{4f'^ + h\Z\\ (1< a, /3 < 2) (112) 



with new constants g\, g2- 

Using the connection between Riesz fractional derivative and its Fourier transform [29] 

(_A^)"/2 ^ (Ki)"/2, (-A||)/^/2 ^ («|)/^/2, (113) 

we obtain from (112) 

= -i^9^ + 9i{-^i.r'^Z + g2{-^\\f'^Z + u:oZ + b\Z\''Z, (114) 

where Z = Z{t, x, y, z). By changing the variables from (t, x, y, z) to {t, ^,y,z), ^ = x — Vgt, and 
using 

(-A,,)/^/^ = JL = ^ (115) 

we obtain from (114) equation 

= gi{-A^r/^Z + 52(-A||)^/2z + uJoZ + b\ZfZ, (116) 

that can be called the fractional nonlinear parabolic equation. For 52 = we get the non- 
stationary FGL equation (fractional NLS equation) suggested in [26]. Let us comment on the 
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physical structure of (116). The first term on the right-hand side is related to wave propagation 
in a media with fractal properties. The fractional derivative can also appear as a result of ray 

chaos [15, 49] or due to a superdiffusive wave propagation (sec also the discussion in [3, 15] and 
corresponding references therein). Other terms on the right-hand-side of Eqs. (114), and (116) 
correspond to wave interaction due to the nonlinear properties of the media. Thus, Eq. (116) 
can describe fractal processes of self-focusing and related issues. 
We may consider one-dimensional simplifications of (116), i.e., 

+ ^^^^^ 

where Z = Z(i, ^), ^ = x — Vgt, or the equation 

dZ d°'Z 

+ ^^^^^ 
where Z = Z{t, z). We can reduce (118) to the case of a propagating wave 

Z = Z{z-Vgt) = Z{ri). (119) 

Then (118) becomes 

d^Z dZ 

for real Z{ri), and c = ivg. This equation takes the form of a fractional generalization of the 
Ginzburg-Landau equation (FGL), when Vg = 0. Eq. (120) differs from the fractional Burgers 
equation [50, 51] in the structure of the nonlinear term. Nevertheless, an analysis similar to 
[50, 51] may be performed to obtain some estimates on the solution [26]. 

5.2 e-expansion for the FGL equation 

Let us consider e-expansion for a particular case of equation (118) when the time dependence 
can be excluded by the replacement 



Z(i,z) = e"'*Z(z). (121) 



It gives the FGL equation 



g-^+OZ + bZ^ = 0, g = gi, Q = LOo + ac (122) 
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for real field Z = Z{z). Again, let us search for a solution in the form 



and 



Z{z) = Zo{z) + eZi{z) + 



d^-'Z d^Z ^2 ^ 



d\z\^-^ dz^ 

similar to section 3.1. In zero approximation Z{z) = Zq{z) which satisfies the equation 

d'^Zo 



dz"^ 



+ uZq + hZ^ = 0. 



As in section 3.3 consider a particular solution of (125): 

Zo{z) = |2a;/6|V2sech(|a;/5|V2 z) . 

In the first approximation 

(PZ 

9-^ + i^bZi + Cu)Zi + gDjZo = 0, 
where Zq = Zq{z) is defined in (126), and 

DlZoiz) = jZl,'\z) + i dyln{y)[Z^'\z - y) - zf\z + y)]. 



The derivatives in (128) are 



2(D 



2uj 



A^jcosh^jAz) - 2) 
cosl? {Az) 



A\mh{Az){<c> - cosh^(yl^)) 
cosh^(y4z) 



(123) 



(124) 



(125) 



(126) 



(127) 



(128) 



(129) 



where A = |a;/g'|^/^. Eq. (127) with respect to Zi{z) has only the integer derivatives and can be 
solved approximately or numerically. The asymptotics z — oo for (122) is discussed in Appendix 
2. 



6 Conclusion 

Equations with fractional derivatives can be considered as a convenient tool to model different 
physical processes with dissipation, persistent memory, and long range interaction. To work 
with such type equations, one needs different approximate schemes such as perturbation theory. 
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asymptotics expansions, etc. In our paper we develop e-expansion for the cases when the order 
of fractional derivatives a can be presented as a = n — e, n = 1, 2, and e <C 1. This expansion 

is not uniform with respect to t ^ 1, and two cases et <C 1, ^ 1, are very different. We 
demonstrate how these expansions work for some examples that have exact a solution (fractional 
linear oscillator) and examples when the e-cxpansion is not trivial (fractional nonlinear oscillator, 
fractional Ginzburg-Landau equation). Fractional derivatives are not defined uniquely, and for 
different types of derivatives the e-expansion is different. 
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Appendix 1: ^-expansion for Riemann-Liouville frac- 
tional derivatives of order a = 2 — e 

Prom the equation 




(130) 



where n — 1 < a < n, follows 




n-l 



■k—a 



(131) 



or for n = 2, a = 2 — £ 



r(-i + £) 



/(o) + 



V(e) 



/W(0). 



(132) 



For £ f ^ 1, we use the expansions 




(133) 




(134) 
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Prom Eq. (132) and (9), we obtain 

oVt'f{t) = f^'\t)+sVlf + ..., (135) 

where 

1^1 fit) = V(0) + r V^'^ (0) + (0) ln(i) + 7/(2) it) + r (r) ln(i - T)dT. (136) 

Jo 

Expressions (135), (136) provides the ^-expansion for Riemann-Liouville fractional derivative. 

Appendix 2: Large x asymptotics for FGL equation 

In section 4, we consider the asymptotic t — oo for Caputo derivative. Let us consider a modified 
equation (118), where we replace the Riesz derivative d"Z/d\z\'^ by the Caputo derivative D'^Z: 

dZ 

i-^ =gD''Z + LOoZ + b\ZfZ, (137) 

where Z = Z{t,z), g = gi, and the boundary condition is z'^\t,Q) = 0. Again, exclude the 
term with first derivative, as in (121), by the replacement 

Z(t,z) = e'''*Z{z). (138) 

It gives modified FGL equation 

gD'^Z + u}Z + b\Z\^Z = Q, (139) 
where a) = cjq + ac, and Z(^)(0) = 0. Substitution of (79) into (139) yields 

h\Z\^Z + ujZ- g ^^^ Z(0) ^ 0. (140) 

Using representation 

Z{z) = R{z)e''f'^^\ (141) 
where R{z) = \Z{z)\ is real and (f){z) = arg(Z), Eq. (140) can be rewritten as 

hR^iz) + QRiz) - ff / ° R(0)e^[^W-'^(^)l « 0. (142) 
r(i — a) 

If the constants b, Co, g are real, then 

bR^{z) + cdR{z) ± ff ^.f " , i?(0) « 0, (P{z) = <^(0) + Trn, (143) 
r(l-a) 
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where n is an integer. The solution of (143) is 



for w 7^ 0, or 



Riiz) ^ ±e . ^ Te{9MR{0)z-^+', {z ^ oo) (144) 

LOL [L — a) 



for o) = 0, 6 7^ 0. Finally, 



^i(i,^) « ±4^-p^2(0,0) T£(5/^)z-2+-e-*Z(0,0), (146) 
a;i (1 — a) 



and 



„l/3 ^— a/3 pjat 

^2(t,^) «^ ±^^:^^—-^R-'/HO)Z{0,0) ^ T£'/'(5/&)'/'^^~'+'^/'i?"'/'(0)e^"*^(0,0). (147) 
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